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Abstract. Data from future high-precision Cosmic Microwave Background (CMB) 
measurements will be sensitive to the primordial Helium abundance Y p . At the 
same time, this parameter can be predicted from Big Bang Nucleosynthesis (BBN) 
as a function of the baryon and radiation densities, as well as a neutrino chemical 
potential. We suggest to use this information to impose a self-consistent BBN prior on 
Y p and determine its impact on parameter inference from simulated planck data. We 
find that this approach can significantly improve bounds on cosmological parameters 
compared to an analysis which treats Y p as a free parameter, if the neutrino chemical 
potential is taken to vanish. We demonstrate that fixing the Helium fraction to an 
arbitrary value can seriously bias parameter estimates. Under the assumption of 
degenerate BBN (i.e., letting the neutrino chemical potential £ vary), the BBN prior's 
constraining power is somewhat weakened, but nevertheless allows us to constrain 
£ with an accuracy that rivals bounds inferred from present data on light element 
abundances. 
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1. Introduction 

Forthcoming experiments on Cosmic Microwave Background (CMB) anisotropies such as 
PLANCKf [1], combined with other astrophysical observations, are expected to provide 
detailed information on the cosmological model which describes the evolution of the 
Universe. Their data will allow us to constrain the parameters of these models with 
an unprecedented accuracy. In fact, to fully extract information from such precise 
experimental data, it is crucial to have detailed theoretical models to compare with, 
possibly using reliable "priors" on cosmological parameters or models which can be 
obtained by independent theoretical tools or experimental data. This helps in reducing 
the effect of parameter degeneracies which typically limits the amount of information 
which can be obtained from data. 

In this paper we present one example of this kind, by considering the impact of 
a detailed estimate of the 4 He mass fraction§ Y p = An Rc /n b on CMB data analyses, 
obtained from an accurate prediction of its value from Big Bang Nucleosynthesis (BBN). 
In this framework, Y p is given as a function of the baryon density and, in more 
exotic scenarios, extra relativistic degrees of freedom and/or non-zero neutrino chemical 
potentials p v . 

The amount of 4 He nuclei produced during BBN plays a relevant role at the epoch 
of recombination, as it is one of the parameters controlling the evolution of the free 
electron fraction [2, 3] (it also affects the later phase of reionisation of the Universe, but 
this effect is extremely small [4]). Thus, the CMB power spectrum which is observed 
today depends significantly on this parameter. Typically, in current CMB analyses (and 
in a number of forecasts) , this parameter is fixed to a reference value Y p = 0.24, suggested 
by independent measurements obtained by studying extragalactic HII regions in blue 
compact galaxies, which are however affected by quite a large systematic uncertainty 
[5,6]. This approach is satisfactory with present CMB data, but it is not fully correct 
as it does not take into account the fact that indeed, Y p is strongly correlated to other 
cosmological parameters, in particular the baryon fraction Q^h 2 and the energy density 
Pn in the form of relativistic species, which is usually parameterised in terms of the 
effective neutrino number N eS = 3 + AiV such that 



with T 7 the photon temperature. The standard scenario of three non-degenerate active 
neutrinos corresponds to AiV = 0.046 due to the effect of non-thermal corrections to 
neutrino distributions during the e + — e~ annihilation phase [7] . 

A different (and more consistent) approach is taken in [8], where Y p is considered as 
a free parameter to be obtained as a result of a likelihood analysis along with the other 

| ESA home page for the PLANCK project: http : //astro . estec . esa.nl/SA-general/Projects/Planck/ 
and Planck-HFI web site: http://www.planck.fr/ 

§ We point out that as usual in the literature we refer to Y p as the " 4 He mass fraction" , though strictly 
speaking it is not correct at the percent level, for it neglects the contribution of nuclear binding energy. 




(1.1) 
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relevant cosmological parameters, or rather fixed up to some given error by adopting a 
gaussian prior. This approach leads to a more correct estimate of the baryon density 
and, as we will also show in the following, avoids biasing the point estimates of some 
cosmological parameters, such as the scalar spectral index n s . 

In the present work we will consider yet a different approach previously suggested 
by Ichikawa and Takahashi [9] , which explicitly takes into account the dependence of Y p 
on VL h h 2 (and possibly also on AN and on the neutrino chemical potentials) as obtained 
from BBN. We use the scale-independent ratio £ = /j, v /kT u , where T„ is the neutrino 
temperature and \i v is the chemical potential, which is assumed to be the same for 
all three active neutrino species, due to the effect of flavour oscillations [10,11]. In 
particular, the Helium mass fraction is considered as a known function Y p BBf * (Tl^h 2 , 
AiV, £) and is thus neither fixed a priori, nor left as a free independent parameter. This 
function can be obtained in the framework of standard BBN theory and is found to be 
a smooth monotonically increasing function of Vt^h 2 and AiV, and decreasing with £. 
In particular, the result adopted in the present paper has been computed by using the 
public BBN code PArthENoPE [12]. This tool provides a careful determination of light 
nuclei abundances, with very small uncertainties. PArthENoPE is the result of a large 
reanalysis, updating the whole nuclear reaction network, and also including the effect 
of radiative corrections to neutron/proton weak processes, and the detailed treatment 
of neutrino decoupling described in [7]. For a comprehensive discussions of these issues, 
see e.g. [13]. The theoretical error on Y p BBN is of the order of 0.0002, thus at the level 
of per mille, over the relevant ranges of (Q^h 2 , AN, £) [13]. The accuracy is ultimately 
limited by the present experimental uncertainty on the neutron lifetime. This tiny error 
essentially amounts to imposing a consistency relation between Y p , Vt\,h 2 , AN and £. 

The goal of this paper is therefore to determine how the assumption of a well 
motivated BBN prior on Y p will affect the estimates of cosmological parameters from 
planck data. Of course, the validity of this method is based on the assumption 
of standard BBN (or degenerate BBN (dBBN) in section 2.3, where £ ^ will be 
considered). It also implicitly assumes that the value of N e g does not change from the 
BBN epoch until last scattering. However, possible future evidence for non-standard 
BBN scenarios or a more involved evolution of relativistic degrees of freedom could 
similarly be accounted for by using a different - yet still theoretically calculable - 
functional form of Y p . 

2. Forecast for planck data 

In our analysis, we compare three different ways of treating the Helium mass fraction: 
• Fixed Y p : 

As in most recent analyses and forecasts of CMB data, Y p is fixed to a value of 
Y p = 0.24. This is also the default value set in CosmoMC/CAMB [14, 15]. Note that if 
we set iVefj to its standard value of 3.046 and assume no neutrino asymmetry, the 
value of 0.24 is not consistent with current bounds on the baryon density. At the 
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moment, CMB data are not very sensitive to Y p , so the bias expected from fixing 
Y p in this way is negligible. 

Free Y p : 

Here, no additional assumptions are made and the Helium mass fraction is kept 
completely free, with a top hat prior ranging from to 1. 

BBN prior on Y p : 

Under the assumption that BBN proceeded in the standard way, and that N e g 
remains constant between BBN and last scattering, we exploit the fact that Y p is 
related to £l h h 2 , AN, and possibly £ if neutrinos have a sizable chemical potential. 
In a first-order approach, one could fix Y p to the BBN prediction Y BBN (£l h h 2 , AiV, £) 
calculated by PArthENoPE. That way, however, one would not take into account the 
theoretical uncertainty in Yp BBN . In order to treat the uncertainty properly, we keep 
Y p a free parameter, but, taking 5Y BBN to be gaussian, add the following term to 
the negative logarithm of the likelihood C of each point in parameter space: 

Since recomputing Y BBN for each point would not be practical, we interpolate 
its value from a pre-computed grid. To account for errors introduced due to 
interpolation, we increase the absolute error on Y p BBN to a{Y p ) = 0.0003. Note 
that we do not employ any additional data here e.g., astrophysical measurements 
of primordial element abundances, that may be subject to large systematic errors. 



2.1. Fiducial data and parameter inference 

Following the method described in detail in [16,17], one can generate a set of mock 
CMB data, using the projected specifications of the PLANCK satellite [1] (see table 1). 
However, for the purpose of forecasting errors, it is sufficient to replace the power 
spectrum of the mock data by that of the fiducial model, which stands for an average 
over many possible mock data sets [16]. The data set comprises the TT- and EE-axiXo- 
correlation spectra as well as the TE'-cross-correlation spectrum for multipoles up to 
t = 2500, and we assume a sky coverage of / s k y = 0.65. In our fiducial model we impose 
the standard BBN consistency relation, spatial flatness, and ignore tensor modes; its 
parameter values are summarized in table 2. 

We then perform the exercise of Bayesian parameter inference for a number of 
models, differing in the number of basic free parameters and the treatment of the 
primordial Helium fraction Y p , as described above. We first work under the standard 
assumption that the chemical potential of neutrinos is negligible (//„ <C kT u ), and 
consider two basic models: 

• A minimal model, with six free parameters (Q^Ii 2 , Q^ m h 2 , H , z re , ln[10 lc Mg], ng), 
inspired by the current "vanilla" model. 
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• An extended model, where in addition to the parameters of the minimal model we 
also vary the neutrino mass fraction f v and the number of extra relativistic degrees 
of freedom AiV. The introduction of /„ is motivated by the observation of neutrino 
oscillations, implying a non-negligible effect of neutrino masses on cosmological 
perturbations [18]. The neutrino fraction and AiV are known to be correlated in 
the analysis of CMB data [19-21]. 

In section 2.3, we will repeat the analysis with one extra free parameter which is 
known to modify the outcome of BBN predictions: a non-zero chemical potential for 
neutrinos. 

We employ a modified version of the Markov- Chain-Monte- Carlo code CosmoMC 
[14] to infer the posterior probability density from the data. Eight Markov chains are 
generated in parallel; their convergence is monitored with the help of the Gelman-Rubin 
i?-statistic [22], and our convergence criterion is R — 1 < 0.02. 

When facing real data, one could be worried that the theoretical prediction for 
the anisotropy spectra might be insufficient. In particular, issues like recombination 
or foreground contamination need to be better understood. Since we use the same 
numerical code for generating and analysing the data, we implicitly assume in this 
forecast that all systematics are perfectly under control. Thus, our inferred parameter 
errors may be slightly optimistic. 

2.2. Standard BBN 

2.2.1. Minimal model The one-dimensional marginalised posterior probabilities for the 
parameters of the minimal model are presented in figure 1. A first striking observation 
is that fixing Y p "incorrectly" to 0.24 leads to a significant bias of up to one standard 
deviation in the point estimates for the baryon density, Hubble parameter, spectral 
index and primordial spectrum normalisation. This is particularly worrisome for ns, 
since bounds on this parameter are often used to constrain inflationary models. We 
therefore strongly recommend not to fix Y p to some arbitrary value, such as 0.24, in any 
analysis of future data. 

The reason for the bias in these parameters are degeneracies with the Helium mass 
fraction; we illustrate these degeneracies in figure 2. Of these four degeneracies, only 

Table 1. List of the experimental parameters assumed for the PLANCK satellite [1]: 
#beam measures the width of the beam, At.p are the sensitivities per pixel and v is 
the center frequency of the three HFI channels least affected by foregrounds. 



iy/GHz 


^beam 


Ax/fjK 


Ap/^K 


100 


9.5' 


6.8 


10.9 


143 


7.1' 


6.0 


11.4 


217 


5.0' 


13.1 


26.7 
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Table 2. In this table we show the free parameters of our model, their fiducial values 



used to generate the data set 


and the prior 


ranges adopted 


in the analysis. 


Parameter 




Fiducial Value 


Prior Range 


UaHS. llldLLL-l LlLTlolLV 


"dm" 


U.ll 


01 — ► QQ 


Baryon density 


n h h 2 


0.022 


0.005 0.1 


Hubble parameter 


h 


0.7 


0.4 -> 1 


Redshift of reionisation 




12 


3 ^ 50 


Normalisation @ k = 0.002 Mpc -1 


ln[10 10 A s ] 


3.264 


2.7 -> 4 


Scalar spectral index 


n s 


0.96 


0.5 -> 1.5 


Helium fraction 


y P 


0.2477 


1 


Neutrino mass fraction 


u 





1 


Number of extra rel. d.o.f. 


AN 


0.046 


-3^4 


Neutrino chemical potential 







-1 -> 1 



Table 3. This table shows the projected absolute errors on the parameters of the two 
models with zero neutrino chemical potential. For all parameters except the neutrino 
fraction we quote the half width of the minimal 68% credible interval [23] , for /„ we 
give the values of the 68% upper limit (the lower limit being zero). The columns 
labelled "free" show the results when leaving the Helium fraction a free parameter, 
while those labelled "BBN" correspond to the results imposing our BBN prior. Our 
results are in very good agreement with those found in [9]. 



Minimal model Extended model 

free BBN free BBN 



n b h 2 


2.2 x 10 


-4 


1.3 x 10 


-4 


2.3 x 10 


-4 


2.2 x 10 


-4 


fldmh 2 


1.4 x 10 


-3 


1.2 x 10 


-3 


4.9 x 10 


-3 


3.1 x 10 


-3 


h 


7.9 x 10 


-3 


5.9 x 10 


-3 


28 x 10- 


3 


24 x 10- 


3 




0.40 




0.39 




0.41 




0.41 




ln[10 10 A s ] 


0.024 




0.015 




0.024 




0.021 




n s 


7.2 x 10 


-3 


3.4 x 10 


-3 


8.2 x 10 


-3 


7.7 x 10 


-3 


Y P 


0.011 




3.1 x 10 


-4 


0.015 




2.1 x 10 


-3 


AN 










0.26 




0.15 




u 










< 0.041 




< 0.035 





the one with Ub = VL h h 2 is "physical", as was already noticed in [8]. The redshift of 
decoupling z*, for instance, crucially depends on the number density of free electrons 
n e . Between helium and hydrogen recombination, this density is related to the baryon 
number density nb(oc c^b) and Helium abundance by n e = n^(l — Y p ). So a shift in 
due to a change in Y p can be reversed by changing the baryon density. The dependence 
of n e on the Helium fraction also plays a role for the diffusion damping scale d, with 
d oc n~ l l 2 . An increase in Y p , for example, means smaller n e and hence a larger d, 
which leads to an additional suppression of power in the CMB temperature anisotropies 
on small scales. Phenomenologically, this signature is similar to tilting the spectrum 
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Figure 1. Marginalised posterior probabilities for the parameters of the minimal 
model. The dashed purple curves correspond to the case where Y p is fixed to an 
"incorrect" value of 0.24, the red curves have Y p as a free parameter and the thick 
black curves correspond to the case with a standard BBN prior. 



of primordial anisotropics (i.e., lowering ns), which explains the degeneracy with the 
spectral index. 

The degeneracies with the Hubble parameter and the normalisation are only indirect 
ones, since these parameters are themselves correlated with the baryon density and the 
spectral index, respectively. Note that the degeneracy with A$ is not invariant under a 
change of the pivot scale. Had we chosen the pivot at a small scale instead of a large 
scale, one would expect a positive correlation instead of an anticorrelation between the 
two parameters. 

A comparison between the run with free Y p and that with a BBN prior shows that 
imposing standard BBN places an extremely strong constraint on the Helium fraction 
if we do not vary AiV. The width of the 68% credible interval on Y p is 6.2 x 10~ 4 , 
i.e., the error is dominated by the theoretical uncertainty in the prediction of Y p BBN 
(see table 3). Essentially, the precise determination of the baryon density from PLANCK 
data will also nail down the Helium fraction. This should not be surprising, given that 
Yp BBN is relatively flat in the direction of u^. Consequently, the expected errors on 
the cosmological parameters hardly differ from the results of an analysis with fixed Y p , 
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n s In [W W A S ] 



Figure 2. Two-dimensional marginalised joint posterior 68%- and 95%-credible 
contours for the minimal model. The thick red contours correspond to results with free 
Y p , the thin black contours represent the result when standard BBN is imposed. This 
plot illustrates the degeneracies of Y p with other cosmological parameters and shows 
how they can be broken by imposing the BBN prior. 

while the errors on the spectral index, the normalisation, the Hubble parameter and 
the baryon density are up to a factor two smaller than in the case with K p as a free 
parameter. 

2.2.2. Extended model Allowing AiV to vary in the extended model slightly weakens 
the constraining power of the standard BBN consistency relation. Since K BBN is 
somewhat steeper in the direction of AiV, the Helium fraction will be allowed to vary 
over a wider range of values. Still, as we can see from figure 3, the error on Y p is still 
significantly improved (the width of the 68% credible interval is 4.1 x 10~ 3 if we impose 
standard BBN, while for a free Y p it is larger by a factor of seven). 

By adding f v and AiV to the free parameters of the model, we naturally introduce 
new correlations which weaken the bounds on the other parameters. The new 
degeneracies (particularly those with AiV) turn out to be more serious than the 
degeneracies with the Helium fraction. In models with free AiV, constraints on the 
other cosmological parameters (apart from = ^W^ 2 and AiV itself) are essentially 
independent on how well the Helium fraction is constrained. Even fixing Y p to the 
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Figure 3. Marginalised posterior probabilities for the parameters of the extended 
model, still with zero neutrino chemical potential. The dashed purple curves 
correspond to the case where Y p is fixed to a value of 0.24, the red curves have Y p 
as a free parameter and the thick black curves correspond to the case with a standard 
BBN prior. 



"wrong" value of 0.24 does not lead to a bias worth mentioning. The only exceptions 
are AiV and the dark matter density (the latter due to a correlation with AiV, see, 
e.g. [23]), for which the error improves between "free Y p " and "BBN prior on Y p " by a 
factor of 1.7 and 1.5, respectively. 

Note that in these runs, not all posteriors peak at their fiducial values. At first 
glance, this may sound strange since we use the fiducial spectra in place of mock data 
spectra. However, the mismatch just reflects the difference between the one-dimensional 
likelihood profile and the marginalised posterior when the posterior is far from gaussian 
in certain directions (like, in our case, the neutrino fraction). 

2.3. Degenerate BBN 

The effect of neutrino chemical potentials on BBN predictions is twofold. It contributes 
to the radiation energy density so the value of iV e fj for three neutrinos with a common 
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value of £ becomes 




) 



i 



) 



(2.2) 



implying a larger expansion rate of the Universe, an earlier weak process freeze out with 
a higher value for the neutron to proton density ratio, and thus a larger value of Y p . 
Subdominant effects on neutrino decoupling [24, 25] and their eventual distribution after 
the e + -e~ annihilation phase [26] are indeed very small and can be neglected. On the 
other hand, a positive £ for the electron neutrino, i.e, a larger number of u e with respect 
to V e , enhances n — > p weak processes compared to the inverse processes, lowering the 
number of neutrons per proton available at the onset of BBN. 

Bounds on the value of the neutrino chemical potential from BBN data analyses 
have been considered by many authors, see e.g. [24-31]. Since flavour oscillations 
enforce the condition of equal chemical potentials for the three neutrino species [10], 
the (common) value of £ is strongly bounded by the neutron-proton beta equilibrium 
and the observed value of 4 He mass fraction. Adopting a conservative error analysis of 
primordial Y p as in [5], one gets —0.04 < £ < 0.07 [31]. This bound can be evaded in 
non-standard scenarios with extra relativistic degrees of freedom contributing to AN, 
which is strongly degenerate with £. 

We now consider the effect of a dBBN prior on Y p in CMB data analysis, in order 
to see what is the impact of the extra parameter £ on future PLANCK estimates of 
cosmological parameters, including of course the value of £ itself. As for the previous 
standard BBN case, we do not use any direct experimental information on Y p . 

2.3.1. Minimal model plus £ We repeat our analysis for the minimal model with one 
extra free parameter £, leaving Y p either free, or imposing a dBBN prior on it. The 
constraint on cosmological parameters is essentially the same in the case "free £ 7^ 
and BBN prior" as in the case "£ = and no BBN prior" . In other words, by assuming 
dBBN instead of standard BBN, the constraining power of the BBN prior on the six 
basic ACDM cosmological parameters disappears. 

However, it is interesting to see how the degenerate BBN prior improves the 
constraint on £ itself. To this end one can compare the results of the "free" and "BBN" 
cases of table 4 (both for the minimal model and the extended one). The projected 
error on Y p does not change noticeably between the two cases, since now the 4 He mass 
fraction is not only determined by the baryon density and AN, but also depends on 
the value of £. Even very small departure from zero allows quite a large variation of Y p , 
therefore there is not much difference in this case between imposing the BBN relation 
and leaving the Helium mass fraction a free parameter. In turn, this strong dependence 
of Y p on £ is also responsible for the significant decrease of the error on £ in the BBN 
case. Indeed, the CMB alone is able to constrain £ only through its contribution to the 
total energy density during radiation domination, see equation (2.2). In this case the 
68% error is er(£) = 0.34, consistent with previous analyses based on the Fisher matrix 
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Table 4. This table shows the same quantities as in table 3 for the degenerate BBN 
scenario (with one extra parameter £ = fj, u /kT u ). Note that /„ = does not lie within 
the 68% credible interval for the extended model with free Y p , so we quote the limits 
of the interval instead of an upper bound. 



Minimal model Extended model 





free 




BBN 




free 


BBN 




n h h 2 


2.1 x 10 


-4 


2.1 x 10" 


-4 


2.2 x 10~ 4 


2.3 x 10 


-4 


fldmh 2 


2.1 x 10 


-3 


1.3 x 10" 


-3 


5.0 x 10~ 3 


4.9 x 10 


-3 


h 


10. x 10 


-3 


7.3 x 10" 


-3 


28 x 10^3 


27 x 10- 


3 


z rc 


0.41 




0.40 




0.43 


0.41 




ln[10 10 A s ] 


0.023 




0.022 




0.023 


0.024 




ns 


7.4 x 10 


-3 


6.9 x 10" 


-3 


8.0 x 10^3 


8.1 x 10 


-3 


Y P 


0.012 




0.010 




0.016 


0.016 






0.34 




0.061 




0.45 


0.093 




AN 










0.27 


0.27 




u 










0.016 -> 0.063 


< 0.039 





approximation [32-34]. When including a BBN prior, the sensitivity of the CMB to Y p 
further reduces the error down to <r(£) = 0.061. Remarkably, this result is comparable 
to the error expected using observations of light element abundances, in case one adopts 
a conservative approach on the Y p error estimate to account for possible systematics. 
We conclude that future CMB data like that from PLANCK will be a very useful probe of 
the neutrino asymmetry. In case of a £ ^ detection, it would be of particular interest 
to compare this finding with constraints from primordial element data. 



2.3.2. Extended model plus £ Lastly, we analyse the extended model with a neutrino 
chemical potential. This requires some particular modifications of CAMB in order to 
explicitly include the chemical potential of neutrinos and anti-neutrinos in the expression 
of the massive neutrino phase-space distribution (as explained in [35]). The total density 
during radiation domination is now parameterised as 

A',,r = 3 [ I — I — I - ] 4(~) ! — • (2-3) 




Our results are shown in table 4. Note that due to slow convergence of the chains for 
these models, we relaxed our convergence criterion slightly, demanding R — 1 < 0.07. 
As for the minimal model, the value of Y p is determined with the same (poor) accuracy 
regardless of whether one imposes the dBBN prior or leaves the Helium mass fraction 
as a free parameter, while the error on £ is strongly reduced in the first case, because 
of the strong dependence of Y p on this parameter. The only other parameter affected 
is the neutrino mass fraction f v . Its degeneracy with £ is illustrated in figure 4. In 
fact, if one were to include large scale structure data the bounds on £ would likely 
be significantly further reduced, because the degeneracy existing between /„ and iV e fj 
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0.1 
0.08 
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t 

Figure 4. Two-dimensional marginalised joint posterior 68%- and 95%-credible 
contours for £ and f u in the extended model with free neutrino chemical potential. 
Red lines correspond to "free Y p " , the thin black lines are the results of imposing the 
BBN prior. 

would be alleviated. Notice also that the error on AN is basically unchanged, as the 
dependence of Y p on this parameter is much weaker. Finally, it is also worth comparing 
the results for the extended models when imposing the standard or degenerate BBN 
priors, i.e. the last column of tables 3 and 4. With the exception of the dark matter 
density, the basic ACDM cosmological parameter errors do not undergo any substantial 
change, while we again see the effect of the introduction of the extra parameter £ in 
reducing the sensitivity to Y p . Correspondingly, the error on AN, and, due to the 
aforementioned correlation, also the error on Q^mh 2 , grow by almost a factor of two 
with respect to the standard BBN scenario, reaching the values which are obtained if 
Yp is taken free parameter and £ = 0, see table 3. 

3. Conclusions 

In this paper we have considered in detail the effects of a careful calculation of the Helium 
mass fraction on future CMB experiments, such as PLANCK, using the known fact that 
Y p is not an independent free parameter, but can rather be fixed in the framework 
of Big Bang Nucleosynthesis as a function of the baryon density, the energy density 
during the relativistic dominated era, as well as, in more exotic scenarios, other physical 
inputs such as neutrino degeneracy. In view of the high precision in parameter estimates 
which is expected to be achieved by PLANCK, this method first adopted in [9] is more 
consistent than the current strategy of fixing by hand the value of Y p . It is also worth 
stressing that it is completely independent from any astrophysical information on light 
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nuclei abundance, as Y p determination from low metallicity HII regions in blue compact 
galaxies. In this respect, this method is different from combined analyses of CMB and 
BBN data (as performed for example in [28,30,36-38]). 

We have considered two different models, the standard BBN scenario, with possibly 
extra relativistic species in addition to three standard active neutrinos, and the case of 
degenerate BBN, with sizable neutrino chemical potentials. 

For standard BBN, Y p only depends on the baryon fraction and AiV. Exploiting 
this functional dependence in a forecast for PLANCK data, we have shown that one can 
avoid a possible bias in the estimate of some cosmological parameters, i.e. the spectral 
index n$ and As, which is instead present if Y p is fixed a priori to some reference value, 
usually given by Y p = 0.24. Furthermore, this method allows for a better determination 
of various parameters (like e.g. the baryon density, the spectral index, the number of 
extra relativistic degrees of freedom and of course Y p itself, see table 3). With a BBN 
prior, the Helium mass fraction can be determined with an accuracy better than 1 %, at 
the level of statistical error of astrophysical determinations, which are however possibly 
plagued by a larger systematic error. On the other hand, without imposing the BBN 
prior, CMB data from PLANCK can determine Y p at the 5-6 % level only. 

In the case of degenerate BBN, due to the strong dependence of Y p on one extra 
parameter - namely, the neutrino chemical potential parameter £ assumed to be flavour 
independent due to flavour oscillations - this result is no longer valid. Imposing the 
BBN prior, the value of Y p is only determined with an order 10 % uncertainty, just as 
if a flat prior was assumed over the whole range < Y p < 1. Nevertheless, exploiting 
the dependence of the Helium mass fraction on £ has a big impact on the way this 
parameter can be determined by CMB anisotropy data. For a fiducial value £ = 0, we 
found that the 68% absolute error on this parameter is 0.06 for the minimal model with 
no extra radiation, and 0.09 for the extended model where AiV and the neutrino mass 
fraction f u are allowed to vary. In both cases, the result is comparable with the error 
obtained when using nuclei abundance data alone. If Y p is assumed to be an independent 
parameter, with no BBN prior, the effect of £ on the CMB power spectrum is only 
via its contribution to the relativistic energy density which shifts the matter-radiation 
equivalence point. Actually, in this case the 68% bound is up to one order of magnitude 
weaker, |£| < 0.45, again for our fiducial model with £ = 0. Further spectroscopic 
measurement of 4 He abundance or a better understanding of systematics effects would 
of course, provide a powerful way of independently constraining (or measuring) the 
lepton asymmetry in the neutrino sector. Yet comparison with future CMB data will 
represent an important consistency check. 
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